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Abstract. We have observed the cluster of galaxies Abell 85 with XMM-Newton. These data have allowed us 
to confirm in a previous paper the existence of the extended 4 Mpc filament detected by the ROSAT PSPC in 
the neighbourhood of this cluster, and to determine an X-ray temperature of about ~ 2 keV. We now present a 
thorough analysis of the properties of the X-ray gas in the cluster itself, including temperature and metallicity 
maps for the entire cluster. These results show that Abell 85 had intense merging activity in the past and is not 
fully relaxed, even in the central region. We have also determined the individual abundances for some iron-group 
metals and a-elements in various regions; the ratios of these metallicities to the iron abundance show that both 
supernova types la and II must be involved in the intra-cluster gas enrichment. Spectral analysis of the central 
region suggests a different redshift of the X-ray emitting gas compared to the mean cluster velocity derived 
from galaxy member redshifts. We discuss the implications of the difference between the cD galaxy redshift, 
the mean galaxy redshift and the hot gas redshift, as well as the possibility of several groups being accreted 
on to Abell 85. Finally, we obtain the dynamical mass profile and baryon fraction taking into account the new 
determined temperature profile. The dynamical mass in Abell 85 has a steep density profile, similar to the ones 
found in A r -body simulations. 



1. Introduction 

Clusters of galaxies are the largest and latest bound, re- 
laxed structures to form in the Universe. They are used 
to probe large scale structure and evolution of galaxies in 
dense environments. A deep understanding of clusters is 
therefore necessary if they are to be effectively used as 
an astrophysical tool and laboratory, and detailed multi- 
' wavelength observations of nearby clusters are one of the 
best means to reach a better knowledge of clusters. 

The complex of clusters Abell 85/87/89 is a very 
well studied struc t ure at both X-r a y [wi th ROSAT, 
IPislar et all (Il997lh iLima Neto et all ljl997lh and with 
BeppoSAX, ILima Neto et alJ l|200lj) ]. and optical wave- 
lengths, with extensive redshift and im aging catalogues 
(Dur ret et alJ ll998a; Slczak et al 1 ll998h ; the mean galax y 
redshift of A bell 85 is z = 0.0555 llDurret et al.lflflQRh^ . 1 
As shown bv lDurret et~aT1 (|l998bj) . galaxies in the region 
of Abell 87 have redshifts comparable to those in Abell 85, 
but Abell 87 is not detected in X-rays (as confirmed by 
the present XMM-Newton data) , suggesting that it is not 
a separate cluster or group. Abell 89 comprises two back- 
ground structures and is not clearly detected in X-rays 



1 At z — 0.0555, 1 arcmin = 90.5h 5 Q kpc, assuming Q,m ~ 
0.3, Ov = 0.7. 



either, so we will not consider these two "Abell clusters" 
any further. 

The general picture drawn from a combined X-ray 
(ROSAT PSPC) and optical (imaging an d large redshift 
catalogue) analysis l|Durret et al.lll998b(l and references 
therein, shows the existence of a main cluster, of a south 
blob which is probably a group at the same redshift as the 
main cluster, and of an extended filament (at least 4 Mpc 
long). The latter may be diffuse emission or a chain of sev- 
eral groups of galaxies, which seem to be falling on to the 
main cluster, the impact region being located in projection 
somewhat north of the south blob. 

In order to analyse better the morphology and physical 
properties of the X-ray gas in various regions of the main 
cluster and in the filament, we have observed Abell 85 with 
XMM-Newton. O ur first results concernin g the filament 
were d escribed by ILima Neto et alJ (|2002^ ; iDurret et alJ 
( 2003) . We present here the results concerning the physical 
properties of the X-ray gas i n the main cluster. They are 
compl ementary to those of iKempner. Sarazin fc Ricked 
(2002), who discuss the Chandra observation of Abell 85 
with emphasis on the X-ray concentration at the north- 
ernmost end of the filamentary structure: the south sub- 
cluster or south blob (the filament lies outside their field). 
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Fig. 1. XMM MOS1 and MOS2 raw image of Abell 85, with the regions analysed in detail superimposed. The concentric 
circles around the cluster centre (white full lines) define the centre and rings used for the detailed spectral analysis; 
the green dashed circle defines the VSSRS region; the yellow dashed ellipse defines the impact region, where the South 
Blob is probably hitting the main body of Abell 85; the blue dotted circle defines the South Blob region. The small 
black circles are point sources that were excluded when accumulating the spectra. The VSSRS region was also excluded 
from the external ring. 
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Fig. 2. Left: Optical "true-color" r-g-i image from the SDSS 2nd Data Release overlayed with logarithmically spaced 
Chandra isocontours. Right: Chandra X-ray image in the [0.3-8.0 keV] band. The small circle on both panels indicates 
an X-ray point source which may be an AGN. The arrow at the bottom right shows the position of a "hole" of X-ray 
emission south of the nucleus clearly visible on the Chandra image. 



Durret et al.: XMM-Newton view of Abell 85 
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2. Observations and data reduction 



Two XMM-Newton l|jansen et al.l200l|) observations were 
performed on January 7th, 2002. The first exposure was 
pointed between the center of Abell 85 (coincident with 
the cD galaxy) and the south blob; the second observa- 
tion was aimed at the southern filament. Both total ex- 
posure times were 12.5 ks using the medium optical filter 
in standard Full Frame mode. The basic data process- 
ing (the "pipeline" removal of bad pixels, electronic noise 
and correction for charge transfer losses) was done with 
SAS V5.3. After the paper was submitted, SAS V6.0 be- 
came available and we tested that several of our results 
(the temperature profile, the South blob spectral fits) re- 
mained within the error bars, thus making a new analysis 
unnecessary. 

For the spectral analysis, we have used the observa- 
tions made with the EPIC/MOS1, MOS2 and PN cam- 
eras. For the EPIC/MOS cameras, after applying the 
standard filtering, keeping only events with FLAG = 
and PATTERN < 12, we have removed the observation 
times with flares using the light-curve in the [10-12 keV] 
band using a 3-sigma clip technique. The cleaned MOS1 
and MOS2 event files have remaining exposure times of 
12 381 s and 12 386 s, respectively. 

We did similarly for the PN, applying the standard 
filtering and keeping only events with FLAG = and 
PATTERN < 4. The cleaned PN event file has a useful 
exposure time of 9 039 s. 

The redistribution and ancillary files (RMF and ARF) 
were created with the SAS tasks rmf gen and arf gen for 
each camera and each region that we analysed. 

The smoothed, merged image of the four cleaned MOS 
observations i n the 0.3-5.0 keV en ergy band has already 
been shown bv lDurret et al.l (20.03), confirming our previ- 
ous detection of the "filament" with the ROSAT PSPC. 

Spectra were analyzed with xspec 11.3. Since the 
spectra were rebinned, we have used s tandard \ 2 mini- 
mizat i on. We have applied both mekal ijKaastra fc Mewel 
1 19931 iLiedahl et alTll995h and APEC (Chandra X-Ray 
Center, ATOMDB Version 1.3.1) plasma models. 

The metal abundances ( or metallicities) are base d 
on the solar values given bv lAnders & Grevessel l|l989l) . 
We chose to use these abundances for easier compari- 
son with previous work, even though these values may 
not be accu rate: more recent solar abundance determi- 
nations (e.g. iGrevesse fc SauvaHl998l IWilms et al.ll2000|) 
give significatively lower O and Fe abundances than 
lAnders fc Grevessel l|l989h . This is important to take into 
account when we use "variable metallicity" models in or- 
der to fit the individual metal abundances (the VMEKAL 
and VAPEC plasma models). 

The photoelectric absorption - mainly due to neu- 
tral hydrogen in our galaxy - was computed using the 
cross- sections given bv lBalucinska-Church fc McCammonl 
(1992), available in XSPEC. 



3. Global properties of the X-ray gas 

3.1. Overall morphology 

The overall XMM-Newton image of Abell 85 is displayed 
in Fig.n with the various regions analyzed in detail below 
indicated. The main cluster, South Blob and VSSR region 
are clearly visible. 

A zoom of the optical SDSS 2nd Data Release 2 and 
Chandra maps of the very central region is displayed in 
Fig. |2 An X-ray point source which may be an AGN is 
visible southeast of the cluster center. The distribution of 
the X-ray gas towards the west and northwest edges of 
the central region seems somewhat more clearcut than in 
the rest of the cluster centre. This suggests a cold front 
comparable to tha t previ ously found in other clusters, e.g., 
iMarkevitch et al ] <l2000l) . A "hole" in X-ray emission is 
also visible south of the cluster center. 

3.2. Temperature and metallicity maps 
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Fig. 3. Temperature map and XMM isocontours in the 
[0.5-8.0] keV band for Abell 85. The scale for the temper- 
ature is keV. 



The temperature and metallicity maps are made in a 
grid, where each pixel is 512 x 512 XMM EPIC physical 
pixels, i.e., each cell grid is 25.6" x 25.6". So, from now on 
"1 pixel" is actually this fat 25.6" x 25.6" "pixel". In each 
pixel we try a spectral fit to determine simultaneously the 
temperature and metallicity. 



Sloan Didital Sky Survey: http://www.sdss.org/dr2/ 
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Durret et aL: XMM-Newton view of Abell 85 
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Fig. 4. Metallicity map with XMM X-ray intensity con- 
tours in the [0.5-8.0] keV band superimposed. The scale 
for the metallicity is in solar units. 

We set a minimum count number (~ 600 counts for a 
temperature map, ~ 4000 for a metallicity map) necessary 
for proceeding with a spectral fit using the MEKAL plasma 
model. 

If we do not have the minimum count number in a 
pixel, we try a square region of 3 x 3 pixels; if we still 
don't have the minimum number of counts, we try a 5 x 5 
pixel region. For the metallicity map only, we went up to 
a 7 x 7 pixel region. If we don't have enough counts, the 
pixel is ignored and we proceed to the next neighbouring 
pixel. This is done for all pixels in the grid, covering most 
of the MOS image. 

When we do have enough counts, the spectral fit is 
done with the hydrogen column density fixed at the local 
galactic value - for each pixel, we estimate JVh using the 
task nh from ftools (wh ich is an interpolation from the 
iDickev fc Lockmanl l|l990j) galactic Nh table). Note that 
we compute the RMF and ARF matrix for each pixel in 
the grid. 

Since the cluster covers more or less the entire XMM- 
Newton field, the background was taken into account by 
extracting spectra (for MOS1, MOS2 and PN) from the 
EPIC blank sky templates described by Lumb et al. (2002) 
and publicly available. We have applied the same filtering 
procedure to the background event files and extracted the 
spectra in all the regions. Finally, the spectra have been 
rebinned with the grppha task, to reach at least 20 counts 
per energy bin. We checked that there was no difference 
in the soft Galactic emission (in the [0.3-1.0] keV range) 
between the source and background files. We also checked 
that the 10-12 keV count rates were the same in our data 
and in the background files. 



The temperature map (Fig. |3J) reveals that the gas is 
colder near the center, as expected from its cooling flow 
classification. Several hotter patches are observed on the 
west half of the cluster. The fact that the region located 
just north of the south blob is hotter agrees with the sce- 
nario in which this is the "impact" region where the groups 
forming the filament hit the cluster, and therefore com- 
press and heat the intracluster gas, probably through a 
shock. 

The very steep spectrum radio s ource (VSSR) 
llBagchi et alJl998tlLima Neto et al1l200l|) located south- 
west of the center corresponds to a somewhat cooler region 
(see its full spectral analysis below). 

Typical errors on temperatures in this map are about 
10%. 

The metallicity map (Fig.fJJ was drawn by setting the 
ratios of the heavy element abundances relative to that 
of iron equal to the solar ratios. Regions in grey are those 
where the spectral signal to noise ratio was not sufficient to 
derive a metallicity. The central region, with an extension 
towards the northeast, has a notably higher metallicity, as 
expected from t he higher density of galaxies which tend to 
increase the heavy element abundances in the intracluster 
medium through supernova winds. A higher metallicity 
is also observed in two zones more or less symmetrically 
located on either side of the central region towards the 
northwest and southeast, at the extremities of the image. 
The northwest region has a cooler temperature, while the 
southeast region has a hotter temperature, so the higher 
metallicities in these two zones cannot be explained by an 
anti-correlation between the metallicity and temperature. 
Typical errors on metallicities in this map are about 10- 
15%. 

4. Radial profiles 

We have determined the radial temperature, metallicity 
and absorbing hydrogen column density profiles using the 
combined MOS1, MOS2 and PN events. We have used 
circular annuli, with a width determined by requiring a 
minimum of about 2000 counts per annulus; this condi- 
tion was obtained by trial and error, requiring a compro- 
mise between spatial resolution and accurate spectral fit- 
ting. All three parameters (temperature, metallicity and 
absorption) were left free to vary. 

4.1. Temperature profile 

The XMM-Newton data points for the gas temperature as 
a function of radius are shown together with the BeppoSax 
and Chandra points in Fig. There is an obvious de- 
crease of the gas temperature towards the center for radii 
below about 2.5 arcmin (250 kpc). The temperature is 
constant between radii of 250 and 450 kpc, then slowly 
decreases outwards. Note that the overall shape of 
the temperature profile remains the same (within 
error bars) whether we exclude the south blob and 
VSSR regions or not. 
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Fig. 5. Gas temperature as a function of radius obtained 
with XMM-Newton (filled circles), Chandra (empty cir- 
cles) and BeppoSax (empty squares). 
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Fig. 6. Gas metallicity profile from o ur XMM-Newton 
data. For comparison we also plot the iDe Grandi et alJ 
(2004) data for "cool-core" and non "cool-core" clusters. 



In order to take the radial variations of the tempera- 
ture into account when deriving the dynamical mass, we 
searched for a mathematical approximation to the tem- 
perature profile. 

As a mathematical approximation to the temperature 
profile, we fit the following empirical law: 



T{R) = T + 2T 



1/2 



l + (R/nY 



(1) 



where ft is a scale parameter and To is the central tem- 
perature. 

The best least-squares fit to the data is obtained for 
fcT = 2.83 ± 0.04 keV and r t = (630 ± 50)^ kpc. Fig. 
shows the best fit temperature profile superposed on the 
data points; the agreement is fairly good. 

We use the profile given by Eq. as an approxi- 
mation to the deprojected temperature profile, T(r). For 
physically realistic gas models, the difference between the 
actual (deprojected) gas temperature and the measured 
emission weighted temperature should be smaller than 
~ 10%, well within the observational uncertainti es (e.g., 
iMarkevitch et alJ Il999t iKomatsu fc Seliakl l200lj) . Using 
Eq. JTJ as an approximation to the true temperature 
profile allows us to compute analytically the dynamical 
mass profile and the corresponding total (dark matter plus 
baryon) density profile (cf. Sect. I6.2|l . 

4.2. Metallicity profile 

The radial metallicity profile is displayed in Fig. [5] The 
metallicity is close to solar in the innermost radial bin, 
and decreases to about O.3Z for radii larger than about 
350 kpc. 

The metallicity profile of Abell 85 has a shape 
similar to the mean metallicity profiles obtained by 
157 Grandi et al ] ll2004h for their sample of "cool-core" 
clusters, but with systematically lower values. The metal 



abundance of Abell 85 is closer to the non "cool-core" 
values. 

4.3. Absorption profile 
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Fig. 7. Hydrogen absorption column profile from our 
XMM-Newton data. The two horizontal dotted lines 
indicate the Galactic absorption columns towards the 
Centre of Abell 85 and towards the S outh Blob obtained 
with t he FTOOLS nh task, based on iDickev fc Lockmanl 
(|l990i see text for det ails); the dot-dashed l ine is the 
galactic Nu value from Frlartmann fc Burton] l|l997l) . cf. 
lArabadiis fc Breemanl l|l999|) . 



The hydrogen absorption column profile is shown in 
Fig. For c omparison we also indicate the galactic value 
derived from lDickev fc Lo ckmanl lll990f) ( using the task nh 
from ftools) and Hartm ann fc Burtonl (I1997T) [using the 
value quoted bv lHartmann fc Burtonl l)l997|) ]. The values 
from Dickey & Lockman are well above ours, while that 
of Hartmann & Burton is in better agreement with ours. 
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Durret et aL: XMM-Newton view of Abell 85 



Though error bars are too large for this effect to be 
clearly significant, we can note that the absorption column 
tends to be smaller than the Galactic value at almost all 
radii, as previously found in a to tally independent a nalysis 
based on ROSAT PSPC data ilPislar et al.lll997l) . Since 
the galactic Nn maps are made on a coarser resolution 
than the scale at which clusters are observed in X-rays, 
the actual Nh in the direction of Abell 85 may be smaller 
than the value quoted in these maps. 

The absorption profile is apparently not flat: there 
may be stronger absorption in the central region, within 
3 arcmin. Such a behaviour was al r eady noticed with 
BeppoSAX data by iLima Neto et all ll200ll) in Abell 85, 
and in other clusters (e.g.. IPointecouteau et alJ 12004 
from XMM data). The origin of such a n absoption 
gradient, if present, is s till controversial l|Allenl |2000: 
IPointecouteau et aLll2004L and references therein); a pos- 
sible hypothesis is absorption due to very cool molecular 
gas and/or dust grains. Notice however that the error bars 
are still too large to confirm the absorption gradient. 



Table 1. Point sources excluded in the various spectral 
fits, sorted by R.A. 

Source number R.A. (J 2000) Dec. (J2000) 
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5. Physical properties and metal abundances in 
various regions 

We now analyze in detail the spectral properties of vari- 
ous regions of the cluster. A list of the X-ray point sources 
excluded in the various fits is given in Tabled all the ex- 
cluded zones around point sources have radii of 12 arcsec, 
except for source 3 which is a particularly bright Seyfert 
(see below). We derive the temperatures and iron abun- 
dances, and wherever possible the abundances of several 
other heavy elements, without fixing their ratios to the 
iron abundance equal to the solar ratios. An obvious re- 
quirement for this is to have enough counts in the spec- 
tra. As an example, there are 63000, 68000, 38000, 12300, 
14800 and 5600 PN counts in the Centre, Ring 1, Ring 2, 
South blob, VSSR and Impact regions respectively (after 
excluding point sources). 

5.1. The Central zone and two concentric rings 

We first analyze three regions drawn in Fig. ^ starting 
with a circular region of radius R = 100 arcsec encompass- 
ing the cluster center (hereafter region "Center"). This 
corresponds more or less to the cooling flow zone, where 
the temperature drops, and therefore a constant temper- 
ature may not be a very good approximation. 

We then extracted two regions corresponding to con- 
centric annuli, where the temperature is expected to be 
roughly constant: one between radii of 100 and 250 arcsec 
(Ring 1), the other between radii of 250 and 400 arcsec 
(Ring 2). Two small circular zones were excluded from 
Ring 2: one of radius 24 arcsec centered on the Seyfert 
galaxy (source 3 in Table 0) and one of radius 80 arcsec 
centered on the VSSR. It was not possible to extract zones 
further out with a sufficient number of counts to estimate 
independently the various heavy element abundances. 



For these three regions, we have simultane- 
ously fit the spectra obtained with all three detec- 
tors (MOS1, M OS2 and PN) with th e VAPEC and 
vmek al models JKaastra fc Mewd Il993t iLiedahl et~ai] 
19951). using the photoelectric absor ption given by 
Bal ucinska-Church fc McCammor] 1 1992ft . For the central 



region, we also fit the spectra with a CEVMKL model, 
which takes into account the existence of a cooling flow. 
The energy ranges taken into account were 0.3-10.0 keV 
for the two MOS spectra and 0.5-8.0 keV for the PN 
spectrum. The Galactic absorption and redshift were left 
free to vary. For several elements, abundances could not 
be estimated, and were frozen to 0.3 times the solar value. 

Fig. [S] shows the spectra for regions Center, Ring 1 
and Ring 2, together with the best vmekal or CEVMKL 
(cooling-flow mekal model) fits and residuals in terms of 
individual % 2 contribution of each energy bin. 

Tabic |2 summarizes the results obtained for the vari- 
ous fits on these three regions. The iron abundance is es- 
timated with a typical precision of 3% in the Center and 
better than 10% in the two Rings. The other elements 
measured in the central zone are: O, Ne, Mg, Si, S, Ca 
and Ni. Mg and S are not detected in the Rings. 

It is interesting to note that the abundances of the 
various metals relative to solar are not at all the same 
from one element to another. In the central region, Fe and 
O are about 0.3 and 0.4Z Q respectively, while Ne and Si 
are somewhat higher and Ni is much higher. 

Clear metallicity gradients are also present, all ele- 
ments but Ca being notably less abundant in the Rings 
than in the Center, suggesting that there is little or no 
radial mixing or convection of the gas. 

Fluxes and luminosities were calculated based on the 
best fit for the three detectors independently. We give in 
Table 0] the averages of these quantities in the [0.5-2.0 
keV] and [2.0-10.0 keV] bands obtained from MOS1 and 
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Fig. 8. Best spectral fits with residues for XMM-Newton 
MOS1 (black), MOS2 (red) and PN (green). Top: central 
region (R < 100") fit with a cooling-flow model (cevmkl); 
middle: Ring 1 fit with a vmekal (100" < R < 250"); 
bottom: Ring 2 fit with a vmekal (250" < R < 400"). 



MOS2. The bolometric X-ray luminosities (integrated in 
the [0.001-100 keV] interval) are also given. Typical errors 
are of the order of 10%. 

Interestingly, the redshift derived from the various 
X-ray gas fits in the central region is between 0.0530 
and 0.0535, corresponding to velocities of cz — 15900- 
16050 km/s. Fig. illustrates the redshift dependency on 
the quality of the fit. It seems clear that the X-ray emit- 
ting gas in the central 4 arcmin is best fit with z — 0.0533. 

On the other hand, the mean optical redshift of 
Abell 85 is 0.0555, correspondin g to a velocity of cz = 
16650 km/s l|Durret et alj ll998ah The mean optical ve- 
locity of the cluster is therefore larger than that of the 
X-ray gas in the innermost 100 arcsec radius region by 
~ 600 km/s. This difference is significant only at a 2a 
level since the MOS (PN) energy calibration is accurate 
to 5 eV (10 eV), translating to a velocity of 250 km/s 
(500 km/s). Since Abell 85 is known to show velocity sub- 
structures, we extra cted from the velocity catalogue of 
iDurret et al. | lll998ah the galaxies in the innermost 100 
arcsec radius region and with velocities in the cluster in- 
terval, as defined bv lDurret et al.1 l)l998bl) . Four galaxies 
correspond to these criteria, and their mean velocity is 
cz =15903 km/s, a value consistent with the X-ray gas 
velocity given above. 
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Fig. 9. Spectral mekal fit within the central 4 arcmin 
using the three XMM cameras, restricted to the [5.0-8.0 
keV] band. Each line corresponds to one of the 3 detectors, 
cf. Fig. [HI Top panel: temperature and metallicity are free 
and the redshift is set to 0.0555. Bottom panel, when the 



redshift is set free, the best fit is z — 0.0533" 



-0.0015 
-0.0009 



(90% 



confidence level). The reduced x is given in each panel. 



Interestingly, the cD appears to have a higher veloc- 
ity cz =16734 km/s [z = 0.0558 ± 0.0002, biuret et all 
l|l998af) ] implying that it is certainly not resting at the 
bottom of the cluster potential well traced by the X-ray 
gas. This velocity difference is in the upper range of cD 
peculiar velocities measured bv lOegerle fc Hilll ( 200l]), in 
agreement with the general for mation scenario f or Abell 85 
involving several mergers fe.g.. lDubinskilll99S|) . 



5.2. The South blob 

The South blob spectrum was obtained using the Chandra 
image to derive the positions of the point sources, which 
are not all clearly visible on the XMM-Newton image, and 
exclude them from the XMM-Newton spectrum extracted 
for the South blob in a 3.0 arcmin radius region. It was 
then fit by MEKAL and apec spectra. The corresponding 
spectrum and VMEKAL fit are shown in the top panel of 
Fig- El and results are given in Tabled The temperature 
is kT=5.75±0.25 keV. Contrary to the regions previously 
considered, where the redshift is left free to vary in the 
fit, the result here is not acceptable, in view of the mean 
optical redshift, which is found to correspond to a velocity 
of 16826 km/s (based on 16 galaxy redshifts). We therefore 
fixed the redshift to the mean cluster value of 0.0555. 



5.3. The VSSR region 

The region called VSSR corresponds to the very stee p 
spectrum radio source MRC 0038-096 ijLarge et alJll99lh . 
Fig. ^2 shows the high resolution 21cm VLA ra- 
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Table 2. Spectral fits of the central region (radius 100 arcsec) with three different models, and of regions Ring 1 
(radius 100 to 250 arcsec) and Ring 2 (radius 250 to 400 arcsec) with two different models. Abundances are in solar 
units and error bars are la. Values of Nh are in units of 10 20 cm -2 . 



Region 


Center 


Center 


Center 


Ring 1 


Ring 1 


Ring 2 


Ring 2 


Model 


vmekal 


cevmkl 


vapec 


vmekal 


vapec 


vmekal 


vapec 


N H 

kT (keV) 
Redshift 

x 2 


Z - JU -0.12 

4 30+ - 04 
u.uoou +0 00Q2 
1741 


9 7n +um 

4- 1 U-0.02 
r 77 +0.15 
°-' '-0.15 

u.uooo_ 0004 
1670 


9 oo+u.uui 
°°— 0.002 
4 ie+0.03 

u.uooo_ 0004 
1755 


1 70+0. 18 

r 7Q+0.09 

°-' u -0.04 
u - uoz ' -0.0009 

1770 


1 qs^+ u -00l 

d - uy -0.08 

u.uooo_ 0006 
1766 


L194± u o: ^ 
6.28l°;i° 

n n5n5+ 0047 

u.uauo_ 0009 
1437 


1-448-0.002 
6.20in! 

0539+ 0004 
u.uaoy_ 0005 

1427 


dof 


1491 


1490 


1491 


1636 


1635 


1346 


1332 


O 

Ne 
Mg 
Si 
S 

Ca 
Fe 
Ni 


n Qr-u.Ub 
u - zo +0.06 

u - oo +0.10 

n 55-°- 06 

27-° 07 

u - z '+0.09 
U.OO +0 24 

0.408;g;Si| 

1 on-0.25 
1 - Oy +0.24 


u - ou -0.05 
n -lo+O.lO 

60+ 06 

n S5+ 08 

n 07+0. 26 
uo1 -0.24 

409+ 012 

i - 04 -0.24 


n oo+O. uu 
U-a8_o.06 
n 70+0. 15 
U - 7 =>-0.10 

0- 17±g;i| 
0.56t°«? 
0.28t°:» 
0.46lg; 2 3 

n zl c:i +0.014 

u.ioi_ 014 

1 9 r+0.19 

1 - ZJ -0.29 


022 +lH189 

u.u^^_ 22 

32+ 012 
u - oz -o.n 

o.i6tS:SS 

o.56±g:ii 
n 971 +0.012 

u - zl ± -0.013 

46+ 27 

u -^ D -0.28 


20 +u n 

u - zu -0.12 

44+ 19 

U - J ' -0.33 

296+ 0012 
u.zyo_ 015 


32 +u ' 18 

U - JZ -0.15 

n 7s+ ' 18 

U-/O_ .l9 

29+ ' 24 

u - z3 -0.24 

70+0-56 
u - ' "-0.50 

236+ 021 

u - zou -0.021 

1 34+0-54 


25 +u ' 24 

1 01+ ' 27 

1 - ul -0.23 

1 0Q+ - 55 

261 + 024 
u ' ZUI -0 .020 

1 07+0.75 

l-J' -0.45 



Table 3. Spectral fits of the South Blob (circle of radius 200 arcsec), VSSR (circle of radius 150 arcsec) and Impact 
regions, with two different models. Abundances are in solar units and error bars are la. Values of Nh are in units of 
10 20 cm- 2 . 



Region 


South Blob 


South Blob 


VSSR 


VSSR 


Impact 


Impact 


Model 


vmekal 


vapec 


mekal 


apec 


mekal 


apec 


N H 

kT (keV) 


042 +u - U4B 

U.U*Z_ 04 2 
r 74 +0.26 
°- ' 4 -0.26 


042 +u - ()44 

u - u ^ z -0.041 
e 7 n+0.25 
u - ' a -0.27 


0.93+^ 
6.26lg; 22 


1 00 +u ' 3U 

1 - uu -0.30 

6 23+ ' 20 

u - zo -0.20 


147+0.548 
"■IK -0.147 
07+O.67 
°-°'-0.65 


0.082in S 8 8 2 
8.841°;" 


Redshift 


0.0555 


0.0555 


0.0555 


0.0555 


0.0555 


0.0555 


x 2 


775 


774 


916 


917 


504 


502 


dof 


758 


758 


882 


882 


533 


533 


O 

Ne 

Mg 

Si 

Fe 


26+' 1 - 24 

u - zo -0.19 
u - ' J -0.30 
u - dz -0.40 

0.25+ ,; 28 

o.i68i8:°s» 


n 4f)+ u ;ib 

u -^ u -0.30 
1 1 1 +0.49 
1-11-0.43 

0.68l°j 2 

U.OO_ 29 

n iqn+ 037 

U. -L»U_q.q3 6 


201 + - 035 

U.ZU1_ 034 


n 990+0.039 

u - iio -0 .038 


220+ 085 

"'""-0.084 


246+ 095 
u.zio_ 094 



Table 4. X-ray fluxes and luminosities in regions Center, Ring 1 and Ring 2. Fluxes are in units of 10 11 erg cm 
s _1 and luminosities in units of 10 44 erg s _1 . 



Region 


Center 


Ring 1 


Ring 2 


South Blob 


VSSR 


Impact 


Flux [0.5-2.0 keV] 


1.35 


1.29 


0.65 


0.19 


0.22 


0.070 


Flux [2.0-10.0 keV] 


1.81 


2.12 


1.14 


0.29 


0.38 


0.145 


L x [0.5-2.0 keV] 


1.88 


1.81 


0.91 


0.25 


0.31 


0.10 


L x [2.0-10.0 keV] 


2.67 


3.10 


1.65 


0.42 


0.54 


0.21 


l>X,bol 


6.17 


6.64 


3.51 


0.92 


1.17 


0.44 



dio isocontours of ISlee et all l(200ll) superposed on the 
MOS1+MOS2 X-ray image. 

This radio source was thought to be associated 
with the sout h-west X-ray excess already observed from 
ROSAT d a ta llLima Neto et alJll997t iBaechi et ai]ll998ft . 
ISlee et all l)200l|h using higher resolution radio data, also 
associated the radio source to the X-ray excess based on 
ROSAT data. With XMM data, however, the relation be- 
tween the X-ray excess and the radio data is not very 
clear, but we can see from Fig. ^] that the high resolu- 
tion VLA radio source does not coincide spatially with the 
X-ray excess. Such a phenomenon could happen if the rel- 



ativi stic electrons blow away some of the X-ra y emitting 
gas l|Boehringer et al.lll993t iFabian et aLll2000jh 

A spectrum was extracted in a circular region of ra- 
dius 2.5 arcmin corresponding to the VSSR (see Fig. 
The best thermal fit is obtained for kT = 6.2 ± 0.20 keV. 
We also tried a double component fit, including a thermal 
component and a power law. The power-law component 
is intended to represent the inverse Compton emission of 
cosmic microwave background photons on the relativistic 
electrons responsible for t he synchrotron rad io emission 
observed in this location (|Bagchi et al.lll998j) . However, 
the non-thermal component is not unambiguously disen- 
tangled from the thermal emission, and the reduced x 2 is 
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0.5 1 2 5 10 

energy [keV] 

Fig. 10. Photon spectra for 3 individual regions. Top: 
the VSSRS region (2.5 arcmin radius region). Middle: the 
elliptical impact region. Bottom: the South Blob region 
(3.0 arcmin radius region). See text for region locations. 



non-thermal X-ray flux is 1.5 x 10~ 12 erg s _1 cm -2 . The 
middle panel of Fig. 1 101 shows the best fit spectrum. 

There are 8 galaxies in the VSSR region with mea- 
sured redshifts, but the corresponding velocities spread 
between 13781 and 17800 km/s with no concentration at 
a specific velocity, so we cannot extract any kinematical 
information. 

5.4. The impact region 

As mentioned above, the region believed to be undergoing 
the impact of the infalling groups on to the main cluster is 
located between the south blob and the cluster center. We 
extracted the MOS1+MOS2+PN spectrum in an ellipse of 
major axis 2.43 arcmin, minor axis 1.33 arcmin and major 
axis position angle aligned with the east-west direction 
(see Fig.EJ. 

The spectrum of this region is shown in the bottom 
panel of Fig.llOland the corresponding results for a MEKAL 
and an apec model are given in Table El As in the two 
previous cases, if the redshift was left free to vary in the 
X-ray spectral fit, we found unreasonably small values, so 
we had to fix the redshift to 0.0555. As expected from the 
temperature map, the temperature in this zone is notably 
higher than in the rest of the cluster: kT = 8.8 ± 0.6 keV, 
in agreement with the idea that the gas in that region 
is compressed and heated by a shock due to the infall of 
groups on to the main cluster. 




Right Ascension (J2000) 

Fig. 11. MOS1+MOS2 X-ray image with the 21cm radio 



isocontours obtained bv lSlee et alj I 200ll) overlayed. The 



radio source seems to be anti-correlated with the west X- 
ray excess emission at a ~ h 41 m 36 s and 5 ~ — 9°22' 
approximately limited by the dashed ellipse. 



5.5. Metallicity gradients 
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Fig. 12. Heavy element abundances relative to iron in re- 
gions Centre (red filled circles), Ring 1 (empty squares), 
Ring 2 (green filled diamonds) and South Blob (empty tri- 
angles). Various models are superimposed: SNII (full line), 
Hypernovae (long dashes), and SNIa WDD1 (dot-dashed), 
WDD2 (dotted) and W7 (spaced dots). 



not significantly better. Assuming that we have detected 
a non-thermal component (with the above caveat), this 



The ratios of the heavy element abundances relative to 
iron in regions Center, Ring 1, Ring 2 and South blob are 
shown in Fig. El Various theoretical metal yield mod- 



10 



Durrct et aL: XMM-Newton view of Abell 85 



els a re superimposed: "Hypernovae" l|Nakamura et alJ 
1200 ll). SN II. and SN la models WDD1, WDD2 and W7 
l|Nomoto etal~lll997t) . 

In none of the four regions does any single model ac- 
count for the heavy element abundances. However, Fig.lT2l 
shows that in the Centre and Ring 1 SNIa seem capable 
of reproducing the observed abundances. This may be ex- 
plained by the presence of a central cD galaxy which, as 
normal ellipticals, is deficient in type II SNe but may have 
a relativ ely high SNIa rate - roughl y 16 type la SN per 
century jCappellaro k TurattdfeoOll) . 

On the other hand, the high Ne, Ca and Ni abun- 
dances measured in Ring 2 require the contribution of 
SNII - this could have happened during the initial burst of 
star formation within the clusters elliptical ga laxies (e.g., 
iMatteucci fc Gibsonlll995t IPeoino et alJl2002jh Note that 
Ring 2 is dominated by the region towards the northwest 
of the cluster which has a notably stronger metallicity. 
Ne and Mg in the South Blob are also closer to those 
due to SNII models while Si is consistent with SNIa mod- 
els. Therefore, it does not seem possible from our data to 
estimate accurately the proportion of SNIa and SNII re- 
sponsible for the enrichment of the intracluster medium , 
as in Ce.g. lDupke fc WhitdEffflOalbl ISanders et alJl2004l) : 
besides being a function of time, this proportion is proba- 
bly a function of pos ition (or local density) in the cluster 
l|Tamura et al.ll2004 . 

6. Mass determination 

6.1. Gas mass 

The gas density distribution is determined from the radial 
X-ray surface brightness profile, Ex, using the fact that: 



nl(r) 



rdr 



R? 



(2) 



where n e is the electron number density of the intracluster 
gas. For an electrically neutral gas, n e rijj and then 
p(r) s» m p n e (r), where m p is the proton mass. 

We hav e used two analytical profiles to describe S(iZ), 
the Sersic dl968l) law: 



E(R) = E exp 



and the /3- model ((Cavaliere fc Fusco-Fernianolll976| ): 

-3/3+1/2 



E(i?) 



1+ ^ 



(3) 



(4) 



where a and R c are the scale parameters, u and /3 the 
shape parameters, and E the intensity normalization. 

The X-ray brightness profile of Abell 85 was obtained 
with the STSDAS/IRAF task ellipse, using the XMM- 
Newton MOS1 exposure corrected [0.3-8.0 keV] band im- 
age. We have masked the CCD gaps and point sources. In 
Fig- El we show the brightness profile together with three 
different fits: a Sersic profile, a /3-model and a /3-model 
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Fig. 13. Emission profile of the X-ray gas obtained from 
the XMM-Newton MOS1 image (the profile from MOS2 is 
litterally indistinguishable from MOS1); three fits are su- 
perimposed: a Sersic law (black), a /3-model (blue dashes) 
and a /3-model excluding the four innermost points (blue 
dot-dashes) . 



excluding the four innermost points, which correspond to 
a known excess emission with respect to a flat core (cf. 
ICerbal et alJlTflfll iLima Neto et nljfmiiil. 

If all the points are included, the Sersic law provides 
a better fit than the /3-model; however, if the four inner- 
most points are excluded, the Sersic and /3-model become 
indistinguishable. 

We have chosen the Sersic profile to describe the emis- 
sivity profile for its simplicit y and capability to d escribe 
well the whole profile (see e.g. lDemarco et al.l2003^ . There 
is an analyti cal approximation to t he deprojection of a 
Sersic profile l|Lima Neto et al. I ll999lh Thus, if the bright- 
ness profile is given by Eq. (0, the electron number den- 
sity is: 



n e (r) = n (r/a) p exp [-{r/a) 1 
2p=l- 0.6097i/ + 0.05563^ 2 , 



(5) 



and, for Abell 85, a least-squares fit of the brightness pro- 
file gives a = 35±7kpc, v = 0.39±0.02 andp = 0.39±0.03. 

In order to estimate the electronic density, no [the 
intensity normalization in Eq. IjBjl]. we integrate the 
bremsstrahlung emissivity along the line-of-sight within 
4 arcmin in the central region and compare the result 
with the flux obtained by spectral fitting of the same re- 
gion (the normalization parameter of the thermal spectral 
model in xspec, which is proportional to v? e ). We thus 
obtain n = (77.0 ± 6.5) x 10~ 3 cm" 3 . 

The gas mass as a function of radius can then simply be 
derived by integrating the gas density in spherical shells: 



M(< r) 



4ira 3 nom p 



p (r\ u 



(G) 
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where 7(2, x) = J e z z t ~ 1 dt is the incomplete gamma 
function. We use p(r) — m p n e (r), supposing the gas is 
electronically neutral and m p is the proton mass. The gas 
mass is shown in Fig. 1141 

At the limiting radius of our data, r = 10 arcmin or 
about 9OO/150 kpc, the gas mass is 3.5 x 1O 13 M . 



R [arcmin] 
6 8 10 12 14 16 




1600 



Fig. 14. Top: gas mass (green dotted line) and dynamical 
mass (black solid line when including the calculation of 
the temperature profile, red dashed line in the isothermal 
case) as a function of radius. The vertical blue dotted line 
indicates the limiting radius of the data. For the isother- 
mal case, we have adopted the mean emission temperature 
kT = 6.5 keV. Bottom: baryon fraction in both cases. 



6.2. Dynamical mass and baryon fraction 

In order to compute the dynamical mass, the tempera- 
ture profile is introduced. Assuming hydrostatic equilib- 
rium and spherical symmetry, the Sersic profile described 
above and the empirical temperature profile, Eq. Qy, the 
dynamical mass is: 



M dyn (< r) = kT r x 3.7 x 10 10 M Q x 
[p+v (x*)"] (1 + x 2 ) (1 + + x 2 ) 



^+3x 5 / 2 (7) 



{1 + x 2 ) 2 

where x = r/r t , a and rt are in kpc and kTg in keV. 

If we assume an isothermal temperature profile, the 
dynamical mass reduces to: 



M d s y ° n = 3.7 x 10 1 



] kTr 



Mr. 



(8) 



with r and a in kpc. 

The dynamical mass is shown in Fig. 1141 as a function 
of radius, both for isothermal and non-isothermal gas tem- 
peratures. For the isothermal case, we have adopted the 



mean emission weighted temperature kT = 6.5 keV. These 
two curves are very similar within the limiting radius of 
our data, but start to diverge further out, showing that 
even when the temperature does not vary very strongly 
with radius (see Fig.EJ, it has a non negligible influence 
on the total mass estimate at large radius. At the last 
measured point, r — 10 arcmin, M dyn = 4.3 x 10 14 M Q 



and Mj y ° n 



4.7 x 10 14 M n 



The dynamical and gas masses allow us to estimate the 
baryon fraction profile in Abell 85. The mass in baryons 
is computed assuming that ~ 16 % is contributed by the 
galaxies, i.e., M ba ryon ~ 1.16M gas ifWhite et alJll993|) . 

The bottom panel in Fig. El shows the baryon frac- 
tion for both isothermal and non-isothermal temperature 
profiles. It is seen to increase strongly with radius within 
the inner 150 kpc, then to remain roughly constant. This 
agrees with previous findings that dark matt er is strongly 
peaked towards the ce nters of clusters (e.g. iGerbal et alJ 
Il992t lAllen et alJl200l|) , in agreement with numerical sim- 
ulations (e.g. Rasia et al. 2004). The baryon fraction at 
r = 10 arcmin is 0.09 and 0.10 for isothermal and non- 
isothermal profiles, respectively. 

6.3. Virial radius and halo density 

Having the dynamical mass we can compute the radius r$ 
that corresponds to 



P(n)/Pc(z) = S ; p(r s ) 



M, 



dyn 



47rr|/3 



(9) 



where p c (z) is the critical density of the Universe at red- 
shift z. For the SCDM and ACDM models, the virial 
radius corresponds to 5 = 180 and 340 re spectively (cf. 
lLacev fc CoIe1ll993UBrvan & Normanlll998h . 

Eq. © cannot be solved analytically when we take 
into account the temperature profile and the dynamical 
mass given by Eq. (0, but it can be solved numerically 
and thus we have: 

8 : 180 200 340 500 2500 
r s [h-i kpc] : 2697 2563 1995 1671 789 

Errors on r$ are about 15% (la confidence level). Taking 
the virial radius as r v i r = 2/i^q 1 Mpc, the radial range 
where we are able to obtain a temperature profile for 
Abell 85 is 0.005 <> r/r vir £ 0.45. 

We also derived the dynamical mass density as a func- 
tion of radius, that can be compared to the density of dark 
matter halos formed in cosmological N-body simulations. 
In Fig. [15Jwe show th e dyna mical mass density compared 
to the iNavarro et all lll997l "NFW" ) "univ ersal" profile 



Moore et all l|l999l "Moore99") 



for dark halos and to the 
profile. 

In the interval 0.1 < r < 10 arcmin, where the temper- 
ature and surface brightness are reliably determined, both 
theoretical profiles agree quite well (the Moore99 profile 
slightly better) with the density profile we obtained for 
Abell 85. 
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R [arcmin] 




Fig. 15. Dynamical mass density as a function of radius 
derived from our data (full line) and compared to a NFW 
model (dashed line) and to a Moore99 profile (dotted line) . 
Pdyn oc r -1 - 9 is a good fit at small radii. The two vertical 
dotted lines indicate the spatial resolution (left) and the 
limit of our data (right). 

The characteristic radius of the NFW profile that best 
fits our density profile is tnfw = 390/i^q 1 kpc; with the 
above determined r v ; r , this implies that the concentra- 
tion parameter is c = r V j r /rNFW ~ 5. This value is in 
agreement with th e cosmological N-body simulations of 
l.ling fc Sutd l|2000h . in a ACDM scenario. 

7. Discussion and conclusions 

Our XMM-Newton data have allowed us to obtain tem- 
perature and metallicity maps for Abell 85 and to derive 
heavy element abundances in various regions. 

The temperatures obtained for the Center and Rings 
(see Table 13 ag r ee ve ry well with those measured by 
Markevitch et alJ l)l998|) from ASCA data in his regions 1 
and 6 -9 respectively. On the other hand, iDonnellv et alJ 
( 2003) find a fairly homogeneous temperature throughout 
Abell 85, except in a small zone southeast of the center, 
and in a large one towards the northwest edge. This dif- 
ference may be due to the lower spatial resolution and 
sensitivity of ASCA. 

The temperatures and/or metallicitics of several re- 
gions seem to indicate that Abell 85 has undergone or is 
still undergoing merging events. However, there is no in- 
tensity enhancement whatsoever in any of the hotter or 
more metal rich regions, suggesting that there is no gas 
density discontinuity. 

We can note that the SE region is the only place where 
kT and Z are not anti-correlated: it is hot and metal rich; 
however this zone is just at the limit where metallicities 
can be measured, so we will not take it into account in the 
coming discussion. 

Numerical simu lations of a clus t er me rger such as those 
displayed, e.g., bv iBourdin et alJ l|2004h in their Fig. 9, 
show that when two clusters merge the gas temperature is 



not affected immediately, but the final temperature map 
can look strikingly similar to the one we have derived, with 
various hotter zones. Different numerical simulations have 
also shown that star formation activity in galaxies was 
only briefly enhanced by a merger, mayb e less than 3 Gyr, 
according to Fig. 4 in Fui ita et al.1 l|l999|) . but then due to 
ram pressure stripping galaxies saw their star formation 
rate (SFR) decrease. These two sets of simulations can 
help us propose a possible scenario that could account for 
what we are seeing in Abell 85. 

First, there appears to be a "filament" extending at 
least 4 Mpc southeast of the cluster an d probably made o f 
groups falling on to the main cluster l|Durret et alJ l2003) 
and hitting it in the Impact region. In favour of this in- 
terpretation are two main facts: i) the temperature of the 
filament is about 2.0 keV, consistent with that of groups; 
ii) the impact region, located more or less half way be- 
tween the cluster center and the south blob, is notably 
hotter (8.8 keV). This is most probably the place where a 
particularly large group hits the main cluster. 

Second, the NW region is cooler and more metal rich 
than the average cluster. A cooler temperature is expected 
in metal-rich regions, but the fact that this zone is more 
metal rich probably implies that it has undergone star for- 
mation with a higher rate than the rest of the cluster, pos- 
sibly triggered by a group merging from the NW towards 
the SE. Note however that the NW zone, as the NE one, 
is just at the limit where metallicities can be measured, 
so we should take these results with caution. 

Third, the NE region has no temperature enhance- 
ment, but has a notably higher metallicity, here also sug- 
gesting metal enrichment due to enhanced star formation 
and/or ram pressure stripping due to a merger. 

Fourth, there are hot zones distributed as an arc from 
south to NW. These could also be the result of a major 
merger coming from the NW, since as mentioned above 
a strong similarity is observed between our X-ray tem- 
perature map an d that of the numeric al simulations of a 
cluster merger bv lBourdin et al] l)2004|) . 

If a cluster has fallen on to Abell 85 from the NW, 
this would be consistent with the existence of a metal rich 
zone in the NW, due to a brief increase of the SFR which 
produced a metal enrichment of the ICM, followed by a 
decrease in the SFR, in agreement with the homogeneity 
of the metallicity map in the rest of the cluster. Such a sce- 
nario suggests that this merger has taken place more than 
4 Gyr ago in order for the SFR to be presently reduced 
(jFuiita et al.lll9 99). If we follow the same line to inter- 
pret the NE zone, where the gas has a high metallicity 
but a temperature similar to the average cluster tempera- 
ture, this would require a second merger coming from the 
NE, but notably more recent (less than 3 Gyr), so that 
star formation is still triggered, thus enhancing the ISM 
metallicity, but no effects are visible yet on the tempera- 
ture. In this scenario, the impact region due to the infall of 
the "filament" mentioned above would then be a relatively 
old merger. 
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Such a scenario appears to be consistent with the var- 
ious heavy element abundances measured in the Center 
and in regions Ring 1 and Ring 2, which imply chemical 
enrichment of the intergalactic gas by both supernovae of 
type la and of type II, and corresponding to different star 
formation rates. 

Therefore, the interpretation of the temperature and 
metallicity maps derived for the X-ray gas in Abe ll 85 
based on numerical s imulations bv lFuiita et alJ 1 1999) and 
iBourdin et al.1 1 20041) leads us to suggest that Abell 85 has 
undergone three mergers, two rather old ones (older than 
4 Gyr) from the NW and South, and a more recent one 
from the NE. Of course, this is only a schematic scenario 
since it neglects galaxy movements within the cluster, and 
assumes that galaxies enrich the ISM at their own loca- 
tion, through enhanced star formation and/or ram pres- 
sure stripping. As pointed out by the referee, we should 
also be aware that repeated AGN outbursts dredging up 
material from the center could also produce comparable 
effects. More detailed modeling obviously requires a better 
knowledge of the galaxy distribution, in particular of emis- 
sion line galaxies. We hope to achieve this shortly through 
multi band imaging allowing to derive spectroscopic red- 
shifts and Ha imaging to draw the distribution of galaxies 
with high star formation rates. 
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